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ABSTRACT 

Context. Gamma-ray bursts (GRBs) being the most luminous among known cosmic objects carry an essential potential for cosmo- 
logical studies if properly used as standard candles. 

Aims. In this paper we test with GRBs the cosmological predictions of the Gurzadyan-Xue (GX) model of dark energy, a novel 
theory that predicts, without any free parameters, the current vacuum fluctuation energy density close to the value inferred from 
the SNIa observations. We also compare the GX results with those predicted by the concordance scenario A-CDM. 
Methods. According to the statistical approach by Schaefer (2007), the use of several empirical relations obtained from GRBs 
observables, after a consistent calibration for a specific model, enables one to probe current cosmological models. Based on this 
recently introduced method, we use the 69 GRBs sample collected by Schaefer (2007); and the most recently released SWIFT 
satellite data (Sakamoto et al. 2007) together with the 41 GRBs sample collected by Rizzuto et al. (2007), which has the more 
firmly determined redshifts. Both data samples span a distance scale up to redshift about 7. 

Results. We show that the GX models are compatible with the Hubble diagram of the Schaefer (2007) 69 GRBs sample. Such ad- 
justment is almost identical to the one for the concordance A-CDM. ^From this particular analysis we can obtain the corresponding 
values of the matter density parameter Q m describing GX models. When the similar procedure is applied to the Rizzuto et al. 
(2007) and Sakamoto et al. (2007) SWIFT satellite data, we verify that the SWIFT sample does not delineate a Hubble diagram 
as clearly as featured by the 69 GRBs sample. 

Conclusions. The analysis of the samples of Schaefer (2007) and those by Rizzuto et al. (2007) and of SWIFT (Sakamoto et al. 
2007) shows that more data and efforts are needed to elucidate both issues: the gamma-ray bursts/standard-candle and lack of a 
theoretical understanding of the physics subjacent to the empirical relations. 

Key words, gamma ray burst, cosmology 



X 



1. Introduction 

Recently Gamma-Ray Bursts (GRBs) were used for 
cosmography aims, for the analyses of the high- 
rcdshift behavior of the A-CDM cosmological scenario, 
as well as for several alternative cosmologies (see 



ISchaefer (20031. 2007); Bloom, et al. 1 (I2003D: Dai. 


et al. 


(2004 


); iGhirlanda, et al. 


(2004); Friedman & Bloom 


(2005 


); iLiang & Zhang ( 


2005. 20061): Xu et al 1 ( 


2005); 


IWang & Dai (20061): Firmani. et al. 1 (20061): and other pa- 



pers). The potential power of GRBs for cosmological stud- 
ies, obviously, resides both in their very high luminosity, 
which the highest among known astrophysical objects, and 
the fact that they undergo practically no extinction over 
cosmological distances. Hence, one can get the possibility 
to trace significantly higher distance scales than it is pos- 
sible via supernovae, that is, to trace very deep in the ex- 
pansion history of the universe. The situation however, not 
as simple. 

Schaefer (2007) has developed a statistical approach 
based on the empirical correlations obtained from several 
observed GRBs characteristics and obtained the Hubble di- 



agram for the GRBs after calibrated them for the concor- 
dance and other c osmological models. This approach has 
been used also in (|Cuesta. et al. I ((2007)). The key issue is 
the use of the empirical relations, e.g. of Ghirlanda, Liang- 
Zhang and others, in the absence of understanding of their 
underlying nature, i. e., the genuine character of the scatter 
in the GRBs luminosity vs. luminosity indicators relations 
and of their mutual links. To this difficulty Schaefer (2007) 
has shown that, although the scatter of each of empirical 
relations can be not as small, their joint action can lead 
to a smaller scatter, useful for probing certain cosmological 
models. Increase of the statistics and deeper studies of the 
systematics and selection effects will certainly increase the 
informativity of Schaefer's approach. 

In the present paper we use the same approach as 
Schaefer (2007), to obtain the gamma-ray burst Hubble di- 
agram for the cosmological models proposed by Gurzadyan 
and Xue (2003). The original motivation for GX cosmologi- 
cal models is the fact that, they predict the current vacuum 
fluctuation energy densi ty close to the valu e inferred from 
the SNIa observations (|Perlmutter I (|1998l ); iRiess. et al. I 
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(fl998h : iPerlmutter I (fl999l) ; iRiess. et al. I (pOOl [2006TH 
without any free parameter. 



2. The Cosmological Models 

Gurzadyan and Xue (2002, 2003) have derived a formula 
for the dark energy, which fits the observed value without 
free parameters 



7T he 1 



7T C 1 



8 L 2 a 2 8 G a 2 : 



(1) 



where ?i is the Planck constant, the Planck length is 

L p = (HG) 1 c~ 3 / 2 , c is the speed of light, and G is 
the gravitational constant. Here a is the upper cutoff 
scale in computation of vacuu m fluctuations and has to 
be clo se to the event horizon (iDjorgoyski fe Gurzadvanl 
(|2006h ). According to (|Zeldovich I ( |l968h ). the vacuum en- 
ergy |T]) corresponds to the cosmological term.GX for- 
mula CD defines a broad set of cosmological models 
(jVereshchagin I (|2006D). For the latter the exist e nce of a 
separatrix was shown (jVereshchagin fc Yegorianl (j2006a|)), 
which divides the space of cosmological solutions into two 
classes: Friedmannian-like with initial singularity and non- 
Friedmannian solutions which begin with nonzero scale fac- 
tor and vanishing matter density. Each solution is charac- 
terized by the single quantity, a density parameter which 
is defined in the same way as in the standard cosmological 
model Cl m = 8 ^Siif" , where [i is the matter density, H is 

J-H 

the Hubble parameter, and index "0" refers to their values 
today. The separatrix is given by 



9. 
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where K = ±1,0 parametrize the spa tial curvature. 
The o rigin of the separatrix was revealed in (|Khachatrvan I 
(2007))), and attributed to existence of invariants in GX 
models. 

Analytical solutions for the GX models both for 
matter density and the scale factor are obtained 
(jVereshchagin fc Yegorian I (|2006b)V It turns out that the 
most simple solutions for the scale factor are again those 
of the separatrix. In one model it is exponential, in the 
others they are polynomials. Vereshchagin and Yegorian 
(2006b, 2008) generalized GX models to include radiation 
and looked for other consequences of the models. 

The predictions of the GX models were shown to 
be compatible to su pernovae and Cosmic M i crowa ve 
Background data in (jDiorgovski fc Gurzadvan I (|2006l )). 
A likelihood analysis of supernovae and radio galaxie s 
data was performed i n (jVereshchagin fc Yegorian I (|2006cf l. 
iKhachatrvan. et al. I (|2007D V 



3. Luminosity Distance Formula for GX-models 

The models are described by the two equations for the 
mass d ensity and scale factor (jVereshchagin fc Yegorianl 
(j2006bh ) 




(3) 



where a dot denotes time derivative. For matter (pressure 
p = 0) we ha ve solution for matte r mass density with GX- 
dark energy (|Khachatrvan I (j2007h ) 



,gx 



na{t)\ c 2 (t) 
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4 J G{t)a 3 (t)' 
where b rn is a GX-invariant for matter. For the scale factor 



aW = c( 'Vw + ^ fc - 



(5) 



The luminosity distan ce cLl is ([Peebles I (|l993); 
iDalv fc Djorgovski I (j2005f n 

d L (z) = a fk(K s ){l + z), 
1 



C(i) dz, 



h(z) = 



H( 



Ho 



1 , a ° 
1 + z = — 



where t\,g , Z £1X6 normalized distance and redshift, respec- 
tively. The function /fc(K s ) is defined as 



{sin(x), k = 1, 
x, k = 0, 
sinh(a;), k = — 1, 



(0) 



here k is the effective curvature K — 7r 2 /3 as in §2$. The 
luminosity distance dL{z) for GX models 



d L (z) = a (l + z)f k y-j= In 
where 



g{z) - 1 5 (0) + 1 



9( z ) = \lp( z + 1 ) + 1 , a 



5 (0) - 1 g[z) + 1 

8TTb? X 



(7) 



3ao 



, [3 = tt z - k . (8) 



For the separatrix a = b^ x — we have a simple equation 
for the luminosity distance 



d L {z) = a (l + z)f k 



hx\z+ 1| 
\/tt 2 — k 



(9) 



4. Gamma Ray Bursts calibrated with GX models 

For the present analysis we benefit of the largest GRBs sam- 
ple having properly determined redshifts and luminosities 
currently available. The first was constructed by Schaefer 
(2007) and includes 69 GRBs whose main observables: 
time lag, variability, peak energy, maximum energy in 7- 
rays and rise time, which were obtained from the GRBs 
data provided by many 7-ray and X-ray satellites, and 
ground-based observatories, as collected in Greiner's home- 
page: http://www.mpe.mpg.de^jcg/grbgen.html The sec- 
ond one is the 41 GRBs sample with firmly determined red- 
shift as collected by Rizzuto et al. (2007) from the SWIFT 
satellite. 



Herman J. Mosquera Cuesta et. al.: Hubble Diagram of Gamma-Rays Bursts 



3 



We performed the calibration procedure (linear regres- 
sion analysis) of five luminosity relations: time lag vs. lu- 
minosity ("Hag — L), variability vs. luminosity (V — L), 
peak energy vs. luminosity (E'pcak — L), peak energy vs. 
geometrically corrected gamma-ray ener gy (-E Dea k — E^ ) 



and risetime vs. lum inosity (trt — L) ijSchaefer I ([2007) ; 
ICuesta. et ai~l (12007ft ). We use the OLS Bisector method 
( Isobe et al. I (|1999T) ) to find a relation between each pair 
of these GRBs observational properties. The best-fit line 
for all luminosity relations is given by the general expres- 
sion log C = a + b logZ, where C is the luminosity and X is 
the luminosity indicator, and a is the intercept and b is the 
slope in each of the calibration plots here presented. Their 
uncertainties are: 



a 



al + (<7 6 Xi) 2 + (ba Xz ) 2 + a. 



2 

sys 



(10) 



where &i og i Ji defines the standard deviation in the lu- 
minosity Li, a a is the standard deviation in the intercept 
a, at, is the standard deviation in the intercept b, a Xi is the 
standard deviation in each Xi variable representing log/, 
and CT sys is the systematic error associated to each lumi- 
nosity Li estimate. The results of the calibration procedure 
are given in the Table [IJ and the plot for all luminosity 
relations is given in the Fig. ([1]). 



Table 1. Calibration Results 



Luminosity Relation 


a 




b 


(Jh 


Osys 


Tlag — L 


52.20 


0.07 


-1.01 


0.09 


0.36 


V-L 


52.41 


0.08 


1.78 


0.19 


0.47 




52.15 


0.05 


1.69 


0.10 


0.41 


-£'peak - E-y 


50.49 


0.05 


1.62 


0.11 


0.21 


r RT -L 


52.45 


0.07 


-1.22 


0.11 


0.47 



with the A-CDM model, for the parameters given in the 
inset. Fig.© shows the Hubble Diagram for the 69 GRBs 
calibrated with the GX models obtained after a slight varia- 
tion of the separatrix solution for each cosmological model, 
i.e., depending on the curvature k = 0,-1. 
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Fig. 3. COLOR-ONLINE Hubble diagram of 69 GRBs cal- 
ibrated with GX models as obtained from the empirical 
relations provided by Fig.-Q] The blue line represents GX 
models calibrated with Q m = 0.66, k = and Hq = 70 
kms _1 Mpc _1 . Orange line is for the concordance cosmol- 
ogy On = 0.27, n A = 0.73. 



Using the same method as Schaefer (2007), we obtained 
the best estimated distance moduli /iJ3 



M = — 

w 



5 

E 

i=i 



1 /of and 



5 

i=l 



(13) 

where the summation is over the relations with available 
data, p,i is the best estimated distance modulus from the 
z-th relation, and <7j is the corresponding uncertainty. Then 
applying the error propagation law to the equation (|13D 
we obtained the standard deviation associated to this best 
estimated as = 1/y/w. 

The Fig.© presents the Hubble Dia gram for the 69 
GRBs calibrated with the GX models and its comparison 



1 For an outlying source of apparent; m, and absolute; M, 
magnitudes, distance estimates are made through the distance- 
modulus: ft = m — M, which is related to the luminosity distance 
cLl (given below in units of Mpc) 



n r° da 
d L = ao(l + x) I — 
J a aa 

through the expression 



M = 51og 10 d L (z)+ 25. 



(11) 



(12) 



5. Results 

5.1. GX compatibility with the 69 GRBs sample 

W e now compare G X models with the data set of 69 GRBs 
of ISchaeferl J2007). We use general least square technique 
for all models with different density and curvature param- 
eters. We find the best fit as Q. m = 2/3, k = 0. The x 2 
value for that parameters is 1.0 37 with 68 DoF. The best 
fit curve is shown in Fig. (4(a)) for the valu e of the cur- 
vature parameter k = 0, while in Fig.( |4(b)"j ) the best fit 
curve is shown for k = — 1. In both figures, which combine 
gamma-ray bursts, supernova type la and radio-galaxies, 
several other HD from the GX model are plotted for dif- 
ferent parameter f2 m , as indicated. Notice that the discrete 
set of data in the curves of Fig.© are due to the lack of an 
extensive study in the parameter space of f2 m . We plan to 
address this point in a forthcoming communication. 

For the same parameters we hav e \ 2 ~ 2 with 264 DoF 
for SN (|Khachatrvan. et al. I (|2007l )L The pair best fit val- 
ues of the matter density parameter are close to the sep- 
aratrix value for k = 0, —1, 1, which once again shows the 
important role of the separatrix in GX models. Indeed, we 
see that the separatrix fits the best the observational data, 
with the shown values of the matter density parameter f2 m 
and the curvature k. 

Concerning the k = — 1 models, let us note that the el- 
lipticity detected in the Cosmic Microwave Background ra- 
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Fig. 5. COLOR-ONLINE X 2 for the Hubble diagrams vs 
the f2 m for k = 0, —1, 1 GX-models. The fit points (circles) 
correspond to the separatrix. 



diation temperature maps is characteristic of photon beam 
motion in hyperbolic spaces (Gurzadyan et al 2005, 2007). 
Thus, it can act as a model-independent indication of the 
non precisely zero curvature of the Universe (Penrose 2005, 
Wiltshire 2007). 

5.2. GX compatibility with the 41 SWIFT Satellite GRBs 
sample 

In this section we use the recently released SWIFT satel- 
lite, BAT instrument GRBs data (Sakamoto et al. 2007), 
namely, a sample of 41 GRBs known to have firmly de- 
termined redshifts (Rizzuto et al. 2007). We performed a 
similar procedure as Schaefer (2007), by constructing only 
three (3) luminosity indicator empirical relations from the 
SWIFT data. In particular by using the two variabilities 
defined by Rizzuto et al. (2007), which were constructed 
taking into account the specific operational characteristics 
of the BAT instrument when analyzing the GRBs mask- 
tagged light curves. Those luminosity vs. luminosity in- 
dicator empirical relations are as follows [see a more de- 
tailed discussion, definitions and references in Rizzuto et 
al. (2007)]: 

— Luminosity vs. Indicator: Vr Variability. — We built 
this relation based on data of Vr variability provided 
by Rizzuto et al. (2007) and SWIFT Flux data from 
Sakamoto et al. (2007). Rizzuto et al. (2007) proved 
that a Poissonian variance describes the statistical fluc- 
tuations of the GRBs light curves (see Apendix A in 
Rizzuto et al. 2007). The main reason is that the count 
rates that were used in the analysis already had their 
background substracted 

— Luminosity vs. Indicator: Peak Luminosity. — This re- 
lation was constructed for each GRBs as given in Table 
1. of Rizzuto et al. (2007), which was obtained after ex- 
tracting the mask-tagged light curve by using a binning 
time of 50 ms in the energy band [15-350] KeV (redshift 
corrected). The 1 s time interval with the highest to- 
tal counts was found, and assumed as the time interval 



corresponding to 1 s peak count rate in the subsequent 
analysis 

— Luminosity vs. Indicator: Vlp Variability. — The same 
procedure as above, but notice that it was supposed that 
no extra-Poissonian variance had to be substracted from 
the already mask-tagged light curves 

Our results for the SWIFT GRBs data analysis are pre- 
sented in the Figure-^ where the three empirical relations 
(6(a) Vlp vs. L, 6(b) Vr, vs. L, 6(c) -Epeak vs. L) described 
above are plotted. Besides, Fig.(6d-6f) present three Hubble 
diagrams obtained by varying the normalization constant 
in the variability relations. Notice that Fig.(6e) combines 
the SWIFT and Schaefer 69 GRBs data. These figures sum- 
marize the whole process of testing the consistency of the 
GX model (for several values of the parameters k and fl m ) 
with SWIFT GRBs data, as shown in Tables EEU The at- 
tentive reader should bear in mind that only three of those 
calibration processes are illustrated in Figure-6. Fig.(6d-6f) 
reveal also the sensitivity of the calibration procedure for 
these definitions of the GRBs variability on the normal- 
ization constant. The combined HD in Fig(6e) appears to 
be still consistent with the GX cosmological model, but it 
is also clear that the SWIFT GRBs data distort the bet- 
ter delineate HD of the Schaefer's 69 GRBs sample. This 
suggests that more data of high redshift GRBs detected 
by SWIFT would be needed to clarify the viability of us- 
ing them as insight into GRBs precision cosmology. In this 
respect, perhaps the use of the SWIFT 77 GRBs sample 
with measured redshifts collected by Butler et al. (2007) 
may be of some help. That sample upon arrival of other 
new datasets will be used in future studies. 



6. Conclusion 

We have used the approach by Schaefer (2007) to obtain 
the gamma-ray bursts Hubble diagram for the cosmological 
GX models with dark energy. We have shown their fit to 
the observational data of GRBs with redshift up to z ~ 6, 
and have obtained the best fit values for the parameter 
il m defining the models, see Table [3J e.g. the best fits for 
GRB69 correspond to tl m = 0.64±0.005, for SN+RG O m = 
0.508 ±0.004. 

Since GX-models predict the current value of dark en- 
ergy defining from SNla data, this fact is remarkable in 
the context of degeneracy between the cosmological param- 
eters. 



Table 2. The confidence levels for the used datasets sep- 
arately and together, with corresponding values of O m ; all 
cases are for k = 0. 



Sample 


x 2 


^2 771 


SN244+RG20 


1.23 


0.508 


GRB69 


1.04 


0.64 


GRB41 


1.41 


0.66 


SN244+RG20+GRB69 


1.58 


0.66 


GRB69+GRB41 


1.16 


0.66 


GRBA11+SN244+RG20 


1.56 


0.66 



The key issue in such an analysis is up to which extent 
the gamma-ray bursts can be used as standard candles, i. e. 
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in revealing of the genuine scatter in the empirical relations 
used for obtaining the Hubble diagram. Given the lack in 
the understanding of the physics of those relations, han- 
dling the new data (Sato et al 2007; Sakamoto et al 2007; 
Rizzuto et al 2007) must need a particular care. For ex- 
ample, recently Campana et al (2007) claimed weakening 
of the Ghirlanda relation (x 2 up to 2 or 3) using a sam- 
ple of 5 bursts observed by SWIFT. Even assuming that 
this result is correct, it will influence our evaluations negli- 
gibly (x 2 = 1-040 instead of 1.037). However, Ghirlanda 
et al (2007) reconsidered the same SWIFT sample and 
showed that, that relation survives practically unmodified. 
Recently also Cabrera et al (2007) have reported the corre- 
spondence of a sample of SWIFT data with those of other 
satellites and, hence, absence of outliers. 

In sum, the preliminary character of testing of cosmo- 
logical models via the gamma-ray burst Hubble diagram 
indeed remains until the understanding of the nature of 
the used empirical relations and availability of more obser- 
vational data. 

We thank the referee for many helpful comments. The 
last three authors are partially supported by INTAS. 

Table 3. Calibration relations for the GX cosmological 
model with parameters k and fl m as indicated. A and B 
represent the linear regression intersept and slope as defined 
in the text, while a a and <jb represent their respective er- 
rors, obtained through the linear regression after using the 
SLOPES method. Here L is the luminosity, i\ a g is the time 
lag, V is the variability, E pea k is the peak (or maximum) 
energy, E 1 is the corresponding energy in gamma-rays, and 
trt is the rise-time. 



Luminosity Relation 


A 


B 




OB 


k = 0, On = 0.6666 










Tlag — L 


52.2 


-1.08 


0.0724 


0.12 


V-L 


52.4 


1.76 


0.102 


0.247 




52.2 


1.71 


0.0548 


0.112 


-Speak" E-y 


50.5 


1.62 


0.0525 


0.119 


trt - L 


52.5 


-1.25 


0.0691 


0.122 



Table 4. Calibration relations for the GX cosmological 
model with parameters k and £l m as indicated. 



Luminosity Relation 


A 


B 


OA 


OB 


k = 0, fi m = 0.633 










Tlag — L 


52.2 


-1.09 


0.0747 


0.123 


V-L 


52.4 


1.71 


0.101 


0.243 




52.1 


1.67 


0.0539 


0.112 


^peak'E-y 


50.4 


1.6 


0.0538 


0.121 


Trt - L 


52.4 


-1.22 


0.0676 


0.12 
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Table 5. Calibration relations for the GX cosmological 
model with parameters k and n m as indicated. 



Luminosity Relation 


A 


B 


OA 


OB 


k = 0, Qm = 0.645 










Tlag — L 


52.2 


-1.09 


0.0733 


0.121 


V-L 


52.5 


1.78 


0.103 


0.248 


Epe&k.-L 


52.2 


1.73 


0.0552 


0.112 


-C/peak - E~y 


50.5 


1.63 


0.0522 


0.119 


trt - L 


52.5 


-1.26 


0.07 


0.123 



Table 6. Calibration relations for the GX cosmological 
model with parameters k and Q m as indicated. 



Luminosity Relation 


A 


B 


OA 


OB 


k = 0, fi m = 0.76 










Tlag ^ L 


52.1 


-1.05 


0.0706 


0.118 


V-L 


52.4 


1.71 


0.101 


0.243 




52.1 


1.67 


0.0523 


0.111 


^peak - E-y 


50.4 


1.6 


0.0538 


0.121 


r RT - L 


52.4 


-1.22 


0.0678 


0.12 



Table 7. Calibration relations for the GX cosmological 
model with parameters k and Q m as indicated. 



Luminosity Relation 


A 


B 


OA 


o B 


k = 0, r^ se p — 2/3 










Tlag — L 


52.2 


-1.08 


0.0723 


0.12 


V-L 


52.4 


1.76 


0.102 


0.247 


Speak-L 


52.2 


1.71 


0.0548 


0.112 


Speak"Ey 


50.5 


1.62 


0.0525 


0.119 


T RT ~ L 


52.5 


-1.25 


0.0691 


0.122 



Table 8. Calibration relations for the GX cosmological 
model with parameters k and Q m as indicated. 



Luminosity Relation 


A 


B 


OA 


OB 


k = -1, « ro = 0.566 










"Hag — L 


52.3 


-1.11 


0.0755 


0.123 


V-L 


52.5 


1.82 


0.104 


0.25 


Ep ea k-L 


52.3 


1.77 


0.0562 


0.112 


-Speak"E~y 


50.6 


1.66 


0.0516 


0.125 


r RT ~ L 


52.6 


-1.29 


0.0724 


0.125 



Table 9. Calibration relations for the GX cosmological 
model with parameters k and Q m as indicated. 



Luminosity Relation 


A 


B 


OA 


OB 


k = -1, Q m = 0.578 










Tlag — L 


52.3 


-1.1 


0.0746 


0.122 


V-L 


52.5 


1.81 


0.103 


0.249 


Epeak-L 


52.2 


1.75 


0.0558 


0.112 


-Speak-Ey 


50.6 


1.65 


0.0521 


0.121 


r RT - L 


52.6 


-1.28 


0.0713 


0.124 



Table 10. Calibration relations for the GX cosmological 
model with parameters k and Q m as indicated. 



Luminosity Relation 


A 


B 


OA 


OB 


k = -1, n m = 0.605 










Tlag — L 


52.3 


-1.09 


0.0732 


0.121 


V-L 


52.5 


1.78 


0.103 


0.248 


E pea ,k-L 


52.2 


1.73 


0.0551 


0.112 


-Speak~E-y 


50.5 


1.63 


0.0522 


0.119 


r RT - L 


52.5 


-1.26 


0.0698 


0.123 
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Table 11. Calibration relations of the SWIFT satellite 
GRBs data for the GX cosmological model with parameters 
k and f2 m as indicated. Several values of the normalization 
constant were used in an attempt to understand the sensi- 
tivity of the resulting HD to such a parameter (see FigJBH, 
(HF)- Also the systematic error is provided. 



Luminosity Relation 


A 


B 


OA 


TB 


Csys 


k = 0, Sl m = 0.666 












V - L(Vr=0.05) 


49.902 


1.98 


0.487 


0.448 


0.729 


V - L(Vr=0.02) 


49.902 


1.98 


0.487 


0.966 


0.966 


V - L(Vlp=0.05) 


52.744 


1.634 


0.222 


0.3 


0.3 


V - L(Vlp=0.005) 


51.11 


1.634 


0.443 


0.352 


0.352 


£ pcak -L(Band=lMeV) 


52.5 


1.961 


0.19 


0.805 


0.805 
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(a) Time lag vs. L (b) Variability vs. L 




Log[Trt(l+z)] 
(e) Rise time vs. L 



Fig. 1. COLOR-ONLINE Results of the calibration procedure. Notice that all relations were corrected to the rest frame 
of the GRB and also by using the luminosity best-fit line obtained from the nonlinear regression method, (a) Time lags 
for 39 GRBs. (b) Variability for 51 GRBs vs. isotropic luminosity, (c) Speak values for 64 GRBs vs. isotropic luminosity, 
(d) -Speak values for 27 GRBs vs. total burst energy in the gamma rays, (e) Rise time for 62 GRBs vs. isotropic luminosity. 
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Fig. 2. COLOR-ONLINE Hubble Diagram of 69 GRBs calibrated with GX models, (a) Case k — and different values 
of Q m . (b) Case k = — 1 and different values of f2 m . The corresponding parameters are indicated in the insets. 




Fig. 4. COLOR-ONLINE a) The Hubble diagram for SNIa, radio-galaxies (RG), and GRBs calibrated both with GX 
model with k = 0, H = 70 kms^Mpc -1 , and A-CDM with H Q = 73 kms~ 1 Mpc~ 1 . SNIa: green points. RG: magenta 
points. GRBS: red points and black error bars, b) Hubble diagram for GX with k = —1 and Hq = 70 kms _1 Mpc _1 
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(c) -Eocak vs. L (d) Hubble diagram 




Redshift - z Redshift - z 

(e) HD for Vr, Vlp with equal normalization constant 0.05 (f) HD for Vr, Vlp with different normalization constant 

Fig. 6. COLOR-ONLINE Results of the calibration procedure for the SWIFT satellite 41 GRBs with redshift firmly 
determined (Rizzuto et al. 2007). Notice that all relations were corrected to the rest frame of the GRB and also by 
using the luminosity best-fit line obtained from the nonlinear regression method, (a) Variability Vlp for 13 GRBs vs. 
isotropic luminosity, (b) Variability Vr for 37 GRBs vs. isotropic luminosity, (c) E pea k values for 41 GRBs vs. peak 
isotropic luminosity, (d) Hubble diagram (HD) from above relations with Vr = Vlp normalization constant 0.02. (e) HD 
combining 41 SWIFT and 69 Schaefer's GRBs samples, using Vr, Vlp variabilities with equal normalization constant 
0.05. (f) Vr, Vlp variabilities with different normalization constant 



